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Abstract 

The complex Kohn variational principle and the (correlated) Hyperspheri- 
cal Harmonics technique are applied to study the N-d scattering above the 
deuteron breakup threshold. The configuration with three outgoing nucle¬ 
ons is explicitly taken into account by solving a set of differential equations 
with outgoing boundary conditions. A convenient procedure is used to obtain 
the correct boundary conditions at values of the hyperradius ~ 100 fm. The 
inclusion of the Coulomb potential is straightforward and does not give addi¬ 
tional difficulties. Numerical results have been obtained for a simple s-wave 
central potential. They are in nice agreement with the benchmarks produced 
by different groups using the Faddeev technique. Comparisons are also done 
with experimental elastic N~d cross section at several energies. 
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One of the main objectives in nuclear physics is the knowledge of the 
nuclear interaction. In practice, the two-nucleon scattering data are used to 
determine the on-shell nucleon-nucleon (NN) interaction. The off-shell prop¬ 
erties of nuclear potentials and many-body force contributions must be tested 
in systems with A > 2. As a consequence, much work has to be devoted to the 
understanding of the three nucleon bound and scattering states. The Faddeev 
theory has been extensively applied to this problem, and, in particular, to the 
study of the n-d scattering under and above the deuteron breakup threshold. 
The corresponding Faddeev equations in momentum space were originally 
solved by Kloet and Tjon [T| for a central s-wave potential. At present, dif¬ 
ferent numerical techniques are available to solve the Faddeev equations in 
configuration and momentum space. In refs. benchmark calculations for 
n-d scattering were given as reference for new techniques. Realistic potentials 
have been used to calculate the n-d scattering cross-section at different ener¬ 
gies For the p-d channel many accurate experimental data are available, 
but the Faddeev approach to this process becomes difficult due to the Coulomb 
repulsion [||-|^. Such a difficulty is not present in the variational technique 
developed in ref. for the N-d scattering below the deuteron breakup. The 
extension of this above the deuteron breakup threshold is the object of the 
present paper. In particular, the asymptotic conditions to describe three out¬ 
going nucleons in the n-n-p and p-p-n states are explicitly taken into account. 

In the following the important aspects of the approach are briefly outlined 
and various results obtained for a simple s-wave central potential are reported. 
More details on the adopted procedures and results for realistic potentials will 
be presented in a forthcoming paper Q. 

Following ref. Q, the wave function of the system is written as a sum of 
two terms 
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The term is solution of the Schroedinger equation in the asymptotic region 
where the incident nucleon and the deuteron are well apart. The Tc term 
must guarantee an accurate description of the system when the three-nucleons 
are close to each other; moreover, for large interparticle separations it has 
to describe the breakup configurations. Tc is written as a sum of channel 
contributions, labelled by the angular-spin-isospin quantum number. The 
associated two-dimensional spatial amplitudes are expanded in terms of the 
Pair Correlated Hyperspherical Harmonic (PHH) basis fl^ . 

The wave function corresponding to an asymptotic state has the 

form 


1=1,3 L US' 

where the summation is extended over the three different choices of the Jacobi 
coordinates (x, y). L is the relative angular momentum between the incident 
nucleon and the deuteron, S is the spin obtained by coupling the spin j = 1 
of the deuteron with the spin 1/2 of the third nucleon and J is the total 
angular momentum of the system. are the collision matrix (5-matrix) 

elements describing the 2 —> 2 elastic scattering. The function 
ingoing (A = in) or the outgoing (A = out) solution of the two-body N- 
d Schroedinger equation in the asymptotic region. These solutions contain 
suitable regularizing factors at small distances. The explicit form of Tc is 


( 2 ) 


^c(xi,yi) = ^(l)a{xi,yi)ya{jk,i) , 
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where the hyperspherical variables p^ = xj + yf and cos (j)i = Xijp have been 
introduced. The pair correlation functions fa are solutions of a two-body 
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Schroedinger-like equation |^] and are included to accelerate the convergence 
of the expansion. At small interparticle distances they take into account the 
correlations introduced by the strong repulsion of the NN potential and go 
smoothly to unity as the interparticle distance increases. The index a. in the 
sum runs over all the channels compatible with the total angular momentum 
J value, and the antisymmetrization and parity conditions. In numerical 
applications the sum is truncated after including all the important channels. 

The quantities to be determined in the wave function specified by eqs.(2-5) 
are the hyperradial functions u'^{p) and the collision matrix elements. They 
will be calculated by means of the Kohn variational principle. 

Below the deuteron breakup the elastic collision matrix is unitary and the 
problem can be also formulated in terms of the reactance matrix, as done in 
Q. Above the deuteron breakup the complex form of the Kohn variational 
principle is better suited to describe three outgoing particles. In this 
case, the 5-matrix and the hyperradial functions are determined through the 
stationary value of the functional 

i < nsj\H - E\^lsj > . (6) 


Let us first point out the main points of the procedure. After performing 
the variation with respect to the hyperradial functions, the following set of 
coupled equations is obtained: 


X! \_^k% ^k’k' + ^k’k' (z’) + ^k’k' ip) '^k' (p) - E>ak{p) ■ 

a',k' ” " 


Details on the explicit form of the coefficients A, B, C, N and the inhomoge¬ 
neous term D can be found in refs. |^TC]. For each asymptotic state 
two different inhomogeneous terms can be constructed in correspondence to 
the asymptotic ^isj functions with A = in or out. Two different sets of 
hyperradial functions are then obtained by solving the system of eqs.(7) for 


( 7 ) 


4 




the two choices of A. In the subsequent step the two sets are combined to 
minimize the functional with respect to variations of the 5-matrix ele¬ 

ments. This is the first order solution; the second order estimate is calculated 
by replacing the first order solution in eq. 

Appropriate asymptotic conditions must be imposed on the hyperradial 
functions u‘^{p) to completely determine the problem. When p —> oo, they 
should vanish below the deuteron breakup, whereas, for positive total energy 
E, they should be proportional to ex.p{i^/Ep). However, this behavior is 
reached only for very large values of the hyperradius |]l^], and, in addition, 
the presence of the Coulomb potential modifies the free outgoing wave also 
at infinity. Therefore, it is convenient to use the asymptotic behavior of the 
coefficients A, B, C, N and D, entering eqs. ( 0 ) , to obtain the solutions 
at large but finite p. Neglecting terms going to zero faster than p~^, the 
asymptotic expression for the set of differential equations can be cast in the 
form 
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The X term originates from the Coulomb potential matrix elements and it 
shows the expected 1/p behaviour. The kinetic energy operator and the nu¬ 
clear potential contribute both to the h term. The final form @ is obtained 
after orthonomalizing the PHH states at p = oo. In the above equation, 
, L = £q, -|- Lq, -|- 2A: -|- 3/2 and Uak{p) are linear combinations 
of the functions The total number of coupled equations in 

eqs. and (HI) is Nf,q, corresponding to all the considered values of a,k. 

For n-d scattering the x term is zero and, if the coupling term h is ne¬ 
glected, the outgoing solutions of (H) are the Hankel functions H^^\Qp). In 
order to take into account the coupling terms, N^q different solutions of eqs. (|8[) 
of the kind 
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where x is the matrix entering eq. (^), are obtained by choosing = 

0) = Saao^kko- The n > 0 coefficients T are determined by recurrence relations 
obtained from eq. as done, for example, in ref. fl^ . 

The solutions of eqs. ( 0 ) are then matched to specific superpositions of the 
functions (p) by imposing the continuity of the logarithmic derivative 

at a given value p = po. The value of the matching radius po is not relevant, 
provided that the asymptotic form (^) is reached, which is rather well verified 
for Po > 80 fm. With such a condition, it has been numerically tested that the 
solutions are insensitive to variation of po; even in the presence of Coulomb 
potential terms. For p —> oo, such solutions evolve as 


Uak{p) - - E , 

^ / a,ao 

qq/co 

corresponding to the correct asymptotic behavior of three outgoing particles 
interacting via long-range Coulomb potentials . From the above equation, 
it results that the Sa^ parameters are just the inelastic 5-matrix elements de¬ 
scribing the 2^3 breakup process. In the n-d case, the hyperradial functions 
asymptotically reduce to Uak{p) 

The NN interaction model considered in the present paper is the s-wave 
potential of Malffiet and Tjon I-III |T^, with the parameter values given in 
ref. Q. Correspondingly, for n-d scattering, the function Tc(x, y) defined 
in eqs. (3-5) includes only £a = 0 channels. Therefore, the total number of 
channels for the states with total spin 5 = ^ (S' = |) is simply two (one). 
On the other hand, the Coulomb potential is active in all the waves, and so, 
in the p-d case, also channels with angular momenta £a larger than zero have 
been included. For central potentials the elastic part of the collision matrix 
does not depend on J and, moreover, S = S', L = L'. Thus, has been 
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expressed in the usual form exp(2i 

The number of hyper spherical states, i.e. of hyperradial functions, in¬ 
cluded in each channel has been increased until the convergence is reached. 
Typically, eight hyperradial functions per channel are enough for a four digit 
accuracy in the phase shift parameters. The rate of convergence for the s-wave 
phase-shifts ‘^60 and ^60 and inelasticity parameters and in function 
of the number Na of considered states (equations) per channel a is shown in 
table 1 at = 14.1 MeV. A similar trend is found at = 42 MeV and 
the corresponding converged results are shown in table 2. For n-d scattering 
a direct comparison with the benchmark results |p can be done, showing a 
very good agreement. The phase shift and inelasticity parameters have been 
also calculated for p-d scattering. Convergence patterns similar to the n-d 
case are obtained, as can be seen in table 1. It can be observed that there is a 
constant difference of about 3° between the p-d and n-d s-wave phase-shifts 
at the energy values considered. This difference is in reasonable agreement 
with the corresponding results obtained by other authors 

To compare with experimental results, the elastic cross section has been 
calculated for several energy values Ej^ of the incident nucleon. Since the 
nuclear potential here considered is only s-wave active, the differential cross 
section starts deviating from the experimental points above E^^ = 20 MeV [||. 
For lower E]\f values, partial waves up to L = 8 give sizeable contributions 
and have been included. The calculated n-d elastic cross section is presented 
in fig.l at neutron energies E^ = 3, 9,18 MeV, together with the experimental 
data of ref. The p-d elastic cross section is given in fig. 2 at proton energies 

Ep = 3,6,9,18 MeV, in conjunction with the high precision data of ref. |17|. 
Despite of its simple form, the MT(I-III) potential reproduces in a reasonable 
way the experimental cross sections. The differences found at small angles 


disappear when more realistic NN interactions are used |]J]. 
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The calculation of the breakup cross section is easily performed once the 
coefficients given in eq. (10) are known. The 5-matrix unitarity imposes 


the following relation between the elastic and inelastic parameters 

ak 

In all the cases we have considered, this relation is well verified numerically 
with a precision of 10“®. An example of n-d and p-d breakup cross sections, for 
the space star configuration, is reported in fig.3. For the sake of comparison, 
the n-d breakup cross section calculated by means of the Faddeev technique 
is reported, as well. From inspection of the fig.3, it can be noticed the good 
agreement between the results of the two n-d calculations. 

The inelasticity parameter goes to one as L increases, as, for ex¬ 

ample, it has been found in refs. |]^Q for n-d scattering. For instance, at 
En = 14.1 MeV, (1 — < 10“^ already for L > 4. A similar behavior 

has been found for the p-d case. The contributions to both elastic and inelas¬ 
tic cross sections of channels with > 0, included in the p-d case, have been 
found to be nearly negligible. 

We conclude with a few remarks. First of all, the Kohn variational prin¬ 
ciple has been successfully applied to treat scattering processes above the 
deuteron breakup threshold. The complex form of the principle is well suited 
to take into account the boundary conditions for the three outgoing nucleons 
and to obtain the second order estimate of the 5-matrix. The expansion of 
the wave function in the PHH basis allows for lowering the number of hyper- 
spherical states to be included. The problem reduces to the solution of a set 
of second-order inhomogenous differential equations with outgoing boundary 
conditions. For p-d scattering the equations are coupled even in the asymp¬ 
totic region due to the Coulomb potential. However, a simple technique can 
be used to calculate the proper boundary conditions for the wave function at 
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values of the hyperradius p ~ 100 fm. The results obtained for the complex 
s-wave n-d phase-shift parameters are in complete agreement with the bench¬ 
mark calculations of ref. |^. The elastic cross sections have been compared 
with the experimental data. For both processes, n-d and p-d, an overall good 
agreement has been observed, with small differences that may be ascribed to 
the rather simple potential used. 

The extension of the method to realistic NN interactions will be the sub¬ 
ject of a subsequent paper. 
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TABLES 


TABLE 1. Real part of the s-wave phase-shift (in degrees) and inelasticity parame¬ 
ter for the doublet and quartet spin states are given as a function of the number of 

hyperspherical functions considered per channel. The incident nucleon energy is = 14.1 MeV. 


Na 





4 

^0 


n-d 

2 

97.96 

0.5093 


67.01 

0.9933 

4 

105.47 

0.4652 


68.88 

0.9788 

6 

105.51 

0.4650 


68.94 

0.9784 

8 

105.50 

0.4649 


68.95 

0.9782 


p-d 

2 

101.16 

0.5430 


70.92 

0.9905 

4 

108.41 

0.4989 


72.53 

0.9801 

6 

108.44 

0.4986 


72.59 

0.9797 

8 

108.43 

0.4985 


72.60 

0.9795 
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TABLE II. Results obtained for the real part of the s-wave phase-shift (in degrees) and 

the inelasticity parameter are given at the specified incident nucleon energies. For the n-d 

process, the benchmarks results of ref. Q are reported for the sake of comparison. 




n-d at En 

= 14.1 MeV 




Vo 



4 

^0 

present 

105.50 

0.4649 


68.95 

0.9782 

Los Alamos 

105.48 

0.4648 


68.95 

0.9782 

Bochum 

105.50 

0.4649 


68.96 

0.9782 



n-d at En 

= 42.0 MeV 


present 

41.33 

0.5026 


37.71 

0.9034 

Los Alamos 

41.34 

0.5024 


37.71 

0.9035 

Bochum 

41.37 

0.5022 


37.71 

0.9033 



p-d at Ep 

= 14.1 MeV 


present 

108.43 

0.4985 


72.60 

0.9795 



p-d at Ep 

= 42.0 MeV 


present 

43.65 

0.5058 


39.94 

0.9047 
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FIGURES 


FIG. 1. Elastic n-d cross section at neutron energies 


are from ref. |^| 


3,9,18 MeV. The experimental data 


FIG. 2. Elastic p-d cross section at proton energies Ep 


data are from ref. |17| 


3,6, 9,18 MeV. The experimental 


FIG. 3. Laboratory breakup cross section for n-d and p-d scattering at Ejv = 14.1 MeV versus 
the arclength S, for the configuration where two neutrons or two protons are detected at angles 
6 i = 51.02, O 2 = 51.02 and ip = 120 degrees. For the sake of comparison with the calculation 
of ref. 1^, only the contributions from s-waves have been included. The solid (dotted) curve 
corresponds to n-d (p-d) scattering. The solid dots are the n-d results taken from ref. Q. 
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